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Abstract 

In this paper, we are interested in a generalised Vlasov equation, which de¬ 
scribes the evolution of the probability density of a particle evolving according 
to a generalised Vlasov dynamic. The achievement of the paper is twofold. 
Firstly, we obtain a quantitative rate of convergence to the stationary solution 
in the Wasserstein metric. Secondly, we provide a many-particle approxima¬ 
tion for the equation and show that the approximate system satisfies the 
propagation of chaos property. 
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1. Introduction 
1.1. The main equation 

This paper is concerned with the long time behaviour and particle ap¬ 
proximation of solutions of the following equation 

dtPt = -p-VqPt+[/3q + A{q) + B * pt{q) - Xz]-'\/pPt+div,,[{Xp+a z)pt]+A^pt. 

( 1 ) 

The spatial domain is with coordinates {q,p,z), with q,p and z each in 
Subscripts as in Vg and are used to indicate that the differential 
operators act only on those variables. The unknown is a time-dependent 
probability measure p: [0, T] —)■ P(R^‘^); A and B are given R'^-to-R'^ maps; 


Preprint submitted to *** 


January I 4 , 2015 




q;,/ 3 and A are given positive constants. Finally, the convolntion B * pt{q) is 
defined by 

B*pt{q)=[ B{q-q')pt{q,p,z)dq'dpdz. (2) 

jR^d 

It is well-known that Eq. ([1]) is the forward Kolmogorov eqnation of the 
following stochastic differential eqnation (SDE) 


r dQt = Pt dt, 

I dPt = Qt dt - A{Qt) dt-B* ptiQt) dt + XZt dt, (3) 

[ dZt = —A Ptdt — a Zt dt + \/2 dWt, 


where {Wt)t>o is the d-dimensional standard Wiener process. 

The law of the R^'^-valned process {Qt, Pt, Zt)t>o evolving according to 
the SDE ([3]) is a solntion of Eq. ([1]) at time t. Systems of the form as in 
(11]) and (|3]) have been stndied by many anthors. When B = 0 and A is a 
gradient of a potential A = W, ([3]) is called a generalized Langevin dynamic 
in (l^ . 1^ . It is also a special case of the so-called Nose-Hoover-Langevin 
dynamic stndied in [l^, 13, BO,!. In these papers (|3]) is considered as a 
stochastic pertnrbation, via the external heat bath variable Zt, of the classical 
deterministic Hamiltonian system for {Qt, Pt)- The generalised Langevin and 
the Nose-Hoove-Langevin dynamics play an important role in molecnlar dy¬ 
namics becanse they maintain canonical sampling of the original Hamiltonian 
system while being ergodic. When B ^ 0, the involntion term B*pt, which is 
often known as Vlasov term or mean-field term in the literatnre, models the 
self-interaction of the nnderlying physical system. In this case, ([3]) belongs 
to the class of weakly interacting particle systems and have been nsed widely 
in applied sciences, especially in plasma physics, see e.g., B H E3 and 
references therein. Combining both cases, the fnll dynamic (E]) can also be 
seen as stochastic pertnrbation (via the angmented variable Z) of Vlasov’s 
dynamic similarly as the generalised Langevin dynamic. This motivates the 
name in the title of the present paper. 

The pnrpose of this paper is twofold. First, we prove the existence and 
nniqneness of a stationary solntion as well as convergence to the stationary 
solntion of other solntions of O- Second, we provide a particle approx¬ 
imation of solntions of ([T]). We address both issnes nsing the stochastic 
interpretation (|3]) of ([I]) and the so-called conpling techniqne introdnced in 
18, 20] and nsed fnrther in B state the assnmptions and main 


resnlts of the present paper. 
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1.2. Assumptions and main results of the paper 

Throughout the paper, we make the following assumption. 

Assumption 1. i? is an odd map on R^'^. Furthermore, A and B are 
Lipschitz, i.e., there are non-negative constants Ca,Cb such that for all 
qi,q2 e 

\Mqi) - Mq 2 )\ < CA\qi - q 2 \, \B{qi) - B{q 2 )\ <CB\qi-q 2 \- (4) 


The assumption that B is odd is motivated by the fact that in applications, 
B is often a gradient of a symmetric interaction potential which depends only 
on the distance of two different particles in the underlying physical system. 
The Lipschitz conditions (jl]) are standard assumptions to ensure the existence 
and uniqueness of the SDE ([3]) with square-integrable initial data, see e.g., 
1^. It also guarantees the well-posedness of Eq. (HJ: it admits a unique 
measure solution pt G V 2 (R.^'^) for any f > 0 provided that the initial data 
Po £ 'P2(R-^'^), where V 2 {B?'^) denotes the subspace of probability measures 
on R^'^ with hnite second moments. 

The hrst result of the present paper is the following theorem concerning 
the existence and uniqueness of a stationary solution as well as convergence 
towards to the stationary solution of all solutions of ([T]). The trend to equi¬ 
librium is measured in terms of the Wasserstein distance IT 2 on V 2 {B?^). 
In particular, W 2 {p,n), with p,n E V 2 {B?^), can be formulated as the in- 
fimum over a set of all coupled random variables that distribute according 
to p and u respectively (more details follow in Definition [12]) . The usage of 
the Wasserstein distance has two advantages. Firstly, it makes use of the 
relationship between the SDE (|3|) and Eq. ([ID, which is the underlying idea 
of the coupling technique, see Section [2J Secondly, the infimum formulation 
of the Wasserstein distance is convenient in estimating upper bounds, which 
is needed to prove convergence, by choosing any admissible couple of random 
variables. 


Theorem 1. Under AssumptionUl for all positive a,j3 and X there exists a 
positive constant po such that, ifO< Ca + Cb < Po, there exist positive 
constant C and C such that 


W 2 ipt,Pt)<C'e^^W 2 ipo,po), for all t>0, (5) 
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for all solutions {pt)t>o o,nd {pt)t>o of ([I]) with respectively initial data po and 
Pq in V 2 (R^'^). Moreover, ([I]) has a unique stationary solution poo and all so¬ 
lutions {pt)t>o converge towards it exponentially in the Wasserstein distance 

W2{pt,Poo) < C'e-^^W2{po,Poo). (6) 


We recall that Ca and Cb are respectively the Lipschitz constants of the 
conhnement A and interaction B forces. The condition that 0 < Ca+Cb < Vo 
in the theorem means that they are linear-like forces, and that the interaction 
is weak enongh. The proof of this theorem is given in Section 13.11 

We now describe the setnp and state the second main resnlt of the 
present paper. Let be independent standard Brownian motions 

on R^, and (Q^, Po\ ^o)*>i be independent random vectors on with 


law Po ^ "/^(R^) and independent of (W*)j>i. Let 


{N) p(iV) y{N). 


following SDE in (R' 


N,N jjl,N 

) M 


r,N,N y1,N 

■ ■ ■ 1 I I 


zf’^)t>o be the solntion of the 


)t>o 


3d\N. 


'dQi'^ = Pf^dt, 

^ dPf^ = -/? Ql’^ dt - dt-j^ E7=i mi'"" - Qi"") dt + X dt, 

' dZf^ = -X Pf^ dt-a Zf^ dt + y/2 dWi, 

(A)i,N pijN yi,N\ _ pj yi\ 

fWo !0 ) ^0 J ~ Ivo)0) 

m 

This is an interacting particle system since there is an interacting term 
jf'12^=1 d^iQt^ ~ Qt^) the SDE for each i. However, the interaction 
is of weak type in the sense that the interacting term for each i depends only 
on the empirical measnre associated to the whole system 


4 W B(Qp - Qp) = [ B(Qp - q)P(dqdpdz), 


where 


N 


p^idqdpdz) := 




i,N^{dqdpdz). 


2 = 1 


( 8 ) 
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On the other hand, we consider the following non-interacting particle 
system, which is a many-particle version of the SDE ([2]) 

'dQi = pidt, 

^ dPl = -jS Ql dt - A{Ql) dt - B * dt + \ Z\ dt, 
dZl = —A PI dt — a Zl dt + \p2 dW^, 

.(Qo)-^0)^o) ~ (Qo)-^0)-^o)) 


where is the distribution of {Ql,Pl,Zl). Since the initial data and the 
driving Brownian motions are independent {Ql,Pl,Zl)t> with i>l are in¬ 
dependent. Furthermore, they are identically distributed and their common 
law evolves according to o. so is the solution of ([T]) at time t with initial 
data po- In comparison with the system in ([7]), the interacting term is re¬ 
placed by i.e., integrating of B{Ql’^ — ■) with respect to instead 

of the empirical measure . 

In the second main result of the present paper, we prove that as N be¬ 
comes larger and larger, the weakly interacting processes {Qt^, Pt^^ Zl'^)t>Q 
in ([7]) behaves more and more like the independent processes {Ql, P]:, Zl)t>o 
in dH]). In other words, the particle system ([7]) satishes the property of prop¬ 
agation of chaos as dehned in 0 , which we recall in Section |2l 

Theorem 2. Let {Qq, Pq, ZQ, i = be N independent H^'^-valued 

random variables with law po G Let {Q]'^, Pj:'^, Zl’^)t>o,i<i<N cmd 

respectively be the solution of (JTj) and (|9]) with initial datum 
(Qq, Pq, Zq), i = 1,..., N. Under AssumptionUl for all positive a, ft, X there 
exists a positive constant rjQ such that, ifQ< Ca + Cb < rjo, then there exists 
a positive constant C, independent of N, such that for all i = 1,... ,N 


sup E 
t>o 


Ql’^ - Ql 


+ 


pi,N _ p. 


+ 




( 10 ) 


As a consequence, the particle system o satisfies the property of propagation 
of chaos. In addition, the law of any {Ql’^, Pt’^, converges to the 

law pt of {Ql, PI, Zl) uniformly in t and in the Wasserstein metric 

SMpW2{pt'^\pt)‘^ <^. ( 11 ) 

t>0 tv 


We stress that this theorem not only proves the convergence but also 
provides a quantitative rate of convergence. This theorem is proved in Section 
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1.3. Organization of the paper 

The rest of the paper is organized as follows. In Section [2l we recall some 
preliminary knowledge. We present the proof of Theorem [1] in Section 13.11 
and that of Theorem [2] in Section 13.21 


2. Preliminaries 

2.1. Wasserstein metric 

We recall that V 2 (R^'^) denotes the space of probability measures fi on 
with finite second moment, i.e., 



[\q\^ + bT + bH Kdqdpdz) < oo. 


The Wasserstein distance W 2 between two probability measures pi,iz ^ V 2 {^^) 
is defined via 


inf ¥.{\Q-Q\^ + \P-P\^ + \Z-zA, (12) 

(Q,P,Z),(Q,P,Z) V y 


where the infimum are taken over all triples (Q, P, Z) and (Q, P, Z) of random 
variables on with laws /i and iz respectively. We note that convergence in 
Wasserstein metric is equivalent to narrow convergence (i.e., tested against 
continuous and bounded functions) plus convergence of the second moments. 
The Wassertein metric plays an important role in many helds of mathematics 
such as optimal transport and dissipative evolution equations. We refer to 
19| and [ij for expositions on the topic. 


More generally, given a positive dehnite quadratic form Q on R^'^, one 
can define a Q-Wasserstein distance W 2 Q{fz,iz) between E P 2 (R^'^) 


via 


W 2 Q{ia,uY= inf m[Q-Q,P-P,Z-Z 

(Q,P,Z),(Q,P,Z} 


(13) 


where the inhmum is also taken over all triples {Q,P,Z) and {Q,P,Z) of 
random variables on R^*^ with laws p and iz respectively. Moreover, since \/Q 
is equivalent to the Euclidean norm \/\q\‘^ + bP + bP) there exist positive 
constants C and C such that 


CW 2 {ti, iz) < W 2 Q(/i, iz) < C'W 2 {fi, iz), (14) 
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for //, G V 2 (R^'^)- 

The underlyii^ idea of the coupling technique introduced in [l8|, [2^ and 


used further in p, is based on flTT|) . First, we hnd a positive dehnite 
quadratic form Q on such that ([3]) is dissipative with respect to Q. Then 
we apply ffTT|) to obtain estimations in terms of the Wasserstein metric. We 
note that the Wasserstein distance can be dehned via a variety of different 
ways, but for the purpose of this paper, the dehnition above will be the most 
useful because of two reasons as already mentioned in the introduction: it is 
facilitated with the coupling technique and is convenient in obtaining upper 
bound estimates. 


2.2. Propagation of chaos 

We now recall the dehnition of propagation of chaos of a particle system 


introduced in 17|. Letting Xl’^ := (Ql’^, , Zf^), and denoting = 

be the join law 




, dx 




dq^’^dp'^’^dz^’^ ior i = 1,N. Let pi \dx 

of ..., X^’^) where for each i = 1,...,X, Xl’^ satishes ([7]). Let 

A: > 1 be a hxed integer, and let (ii ,... ,ik) be any A;-uple of [1,X]. The 
(* 1 ,..., 4)-niarginal Qf pX) jg dehned via 

p(*i’■■■’*'=) (A) := pI^\A X 


for any Borel set A C We recall that pt is the law of {Qj, Pf, Zj), 

which solves Eq. ([T]). The particle system ([7]) is said to satisfy the property 
of propagation of chaos with respect to pt if converges narrowly to 

{pt)®^ for all fc-uple (ii,... ,4)- According to [13, Proposition 2.2], this is 
equivalent to the statement that the empirical measure p^ converges nar¬ 
rowly to pt- We refer to 0 for more information on the topic. 


3. Proofs of the main theorems 


The proofs of Theorem [T] and Theorem [2] are presented re^ectively in 
Section IXn and Section 1X21 The proofs follow the procedure of |y,l^ (see also 
20| ) which use the stochastic interpretation (E]) of ([I]) and the coupling 


technique explained in the previous section. 


3.1. Stationary solution and long time behaviour of the main equation: proof 
of TheoremUl 

The main ingredient of the proof of Theorem [1] is the following proposi¬ 
tion, which is the content of the coupling technique. 
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Proposition 1. Under the assumption of Theorem [T], there exists a positive 
constant C and a positive dehnite quadratic form Q on such that 


W2Q{pt,Pt) < e ^^W2 q{po,Po), for all t > 0 . ( 15 ) 

Proof. Let {Qt, Vt, Zt) and (Qt, Pt, Zt) be two -valued stochastic processes 
evolving according to ([ 3 ]) with the same Brownian motion {Wt)t>o in Let 
{qt,Pt, Zt) := {Qt — Qt,Pt — Pt-,Zt — Zt) be the difference between them. Then 
{qt,Pt, Zt) evolves according to 

dqt = Ptdt, 

dpt = -jd qt dt - {A{Qt) - A{Qfj dt - {B * pt{Qt) - B* PtiQt)) dt+ \zt dt, 
dzt = —Xpt dt — a Zt dt. 

(16) 

Let Oi, • • • , 05 be positive constants, which will be specihed later on. Dehne 
Q{q,p, z) = ai|g|^ a2|p|^ 204^ ■ p + 205^ ■ z + 2 p- z. ( 17 ) 

Assume that {qt,Pt,Zt) evolves according to flTbl) . Then we have 
d 

^Q{qt,Pt,Zt) 

= 4(ai|gtp a2\pt\^ + a^lztl"^ + 2aiqfPt + 2a^qt ■ Zt + 2pt ■ Zt) 

dt 

= 2aiqt ■ pt + 2a2Pt ■ [-/3 qt - {A{Qt) - A{Qf) - {B * pt{Qt) - B^ 'pti.Qt)) + Zt] 
+ 2a3Zt ■ [-Xpt - azt] + 2a4|pip 

+ 204 ?* • [-/d qt - {A{Qt) - A{Qf) - {B * pt{Qt) - B* pt{Qt)) + X Zt] 

+ 2o5P* • Zt + 2o5?* • [-Xpt - azt] + 2 pt ■ [-Xpt - azt] 

+ 2zi • [—/d qt — {A{Qt) — A{Qf) — {B * Pt{Qt) — B * PtiQt)) + 

= ( 2 ai - 2A05 - 2 a 2 (d)qt ■ Pt + (2A04 - 2ao5 - 2 (d)qt ■ Zt 

P (2A02 — 2A03 -|- 205 — 2 a) Pt ■ Zt — 2/3o 4 |?*P — (2A — 2o4)|p*p 

— ( 2 q ;03 — 2A)|zt|^ — ( 2 o 2 P* -|- 2 o 4 ?* -t- 2 zt)iAiQt) — AiQf)) 

- { 2 a 2 Pt + 204 ?* + 22 :*) (R * PtiQt) - B* PtiQt))- ( 18 ) 


Now we estimate the last two terms on the right-hand side of flTSl) using the 


assumptions on A and B. The first term is bounded from above by 

— {2a2Pt + ^a^qt + 2zt){A{Qt) — A{Q^)) 

< 2a2CA\qt\\pt\ + 2a4^Cjilqtl"^ + 2C'^|gt| 

< a2CAi\qt\‘^ + Iptl"^) + 2a4CA\qt\‘^ + 

= (a2 + 204 + ^)CA\qt\‘^ + ci2CA\Pt\‘^ + (19) 

The second term is a bit more intricate. Let 11^ be the joint law of {Qt, Pt, Zp, Q^, Pt, Zt) 
on X Then its marginals on R^*^ are distributions pt and of 

{Qt,Pt,Zt) and {Qf-, Pt, Zt) respectively. We have 
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Similarly, we have 


- 2E [pt ■ {B * pt{Qt) - B* ptiQt))] 

=-2 [ {P-P).{B{Q-q)-B(Q-q)) dUt{q, p, z; q, p, z) dUt{Q, P, Z; Q, P, Z) 

jRl2d 

= - [ {iP-p)-(P-p))-{B{Q-q)-B(Q-q)) dU,{q, p, z; q, p, z) dn,{Q, P, Z; Q, P, Z) 

J-Rl2d 

<Cb [ \{P - P) - (P -p)\ ■ \{Q - q) - (Q -q)\ dUtiq,p, z; q,p, z) dUt{Q, P, Z; Q, P, Z) 

jRl2d 

< % / [\iP-p)-(P-p)\" + \iQ-q)-(Q-q)\"] dMq, p,p q, p,z) du^iQ,p,Z; g, p,z) 
= % / [\{P-P)-{p-p)\^ + \{Q-Q)-{q-q)\^]dnM,P. q,P, z) dU.iQ, P, Z; Q, P, Z) 

Z Jlll2d 

<CBE{\q,\^ + \p,\^), 
i.e., 


- 2E [p, • (p * p,{Q,) -B*, fiMQt))] < CBH\qt\" + btH- (21) 

In the same way, we also obtain 


- 2E • (P * pM) - B * p,(g,))] < CBE{\q,\^ + \z,\^). (22) 

Substituting estimates from flT^ to fl 2 ^ into flT^ . we get 


d 

■^Q{qt,Pt, Zt) < (2ai - 2Aa5 - 2a2^) qfPt + [2Aa4 - 2aa5 - 2/3] qt ■ Zt 

+ [ 2 Aa 2 — 2 Aa 3 + 205 — 2cx\pt ■ Zt — [2/3a4 — (02 + 204 + V){Ca + C*r)] l^tl^ 
— [2A — 204 — a2(C*A + Cb )) \ pt \^ ~ [20:03 — 2A — { Ca + C*r)] \ zt \^■ 

(23) 


Set p := Ca + Cb- We choose oi,..., 05 such that 


"2ai — 2Aa5 — 202 /? = 0, 

2Aa4 — 2005 — 2/5 = 0, 

2X0,2 — 2Aci3 -|- 205 — 2o = 0, 
2/5ci4 — (02 P 2ci4 P l)h ^ 0, 
2A — 204 — 02 ^ > 0, 

^2003 — 2A — r/ > 0, 


( 24 ) 
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and that Qiq,p,z) is a positive definite qnadratic form on Firstly, we 
choose 04 = I and express Oi, 02 , 05 in terms of 03 . Then we show that there 
exists 03 and r]o > 0 snch that all the reqnirements are fnlfilled for 0 < rj < r]o. 
Conditions in fl 2 T|) are eqnivalent to 

a a ’ 

^2 = 0^3 + f ~ 

Oi = Aos + j3a2 = Aos + /? (^03 + y ~ j 
< ^ ^ (25) 

2/^04 — (ci2 T 204 1 ) T ] = X (3 — (03 + y — 

2A — 204 — 02 ?^ = A — ^03 + f — > 0, 

2aa3 — 2A — 77 > 0. 

Let 03 = 2 + 03 , where 03 will be chosen later on. Then we have 

Q{q,P,z) = (a5{X - j) + ^{as + + (®3 + y - y) |pp + Oskl^ + Ag ■ p 

— 205 q ■ z + 2p ■ z 

= ^a5(A “ y) + /^(®3 + ^ + ^)) + (®3 ^ + (2 + a3)\z\ 

+ Xq ■ p — 205 q ■ z + 2p ■ z 

= (^ 05 (A “ y) + ^(®3 + ^ + 2) - Og - y j \q\‘^ A ~ t) 

/A \ 2 

+ [^q + P) + {a5q +zf + {p +zf. 

Now we choose 03 snch that 

■= o,5{^ “ a) ■*" <^(“3 + f + 2 ) — oi — X > 0 ) 

/ 2 (h 3 ) := 03 + f “ ^ > 0, 

fsi^s) ■= ^3 > 0 , (26) 

74 ( 03 ) :=a 3 + y + A + 3 — ^> 0 , 

^ 75 (^ 3 ) •= 2a:(2 + 03 ) — 2A > 0. 

Since all 7n with z = 1,..., 5, are linear fnnctions of 03 with positive slopes, 
it is obvions that there exists 03 that is large enongh so that all of them are 
positive. We choose snch a 03 and define 

{yM'2 + /2(a3)"^*<“=‘'} 


Po := mm 


(27) 











It then follows from (12^ and (HTh that fl2S]) are fulhlled for every 0 < r] < tjq 
and Q is a positive dehnite quadratic on 

From fl23|l it follows that there exists a positive constant C, depending 
only on a, I3, A, Ca and Cb such that 

^EQ{qt,pt,Zt) <-CE[\qt\‘^ + for all f > 0. (28) 

Since in a hnite dimensional space, all norms are equivalent to the Euclidean 
norm, the right hand side of (12511 is bounded above by —CE Q{qt,Pt, Zt) for 
some positive constant C, i.e., 

^EQ{qt,pt,Zt) <-CEQ{qt,pt,Zt), for all f > 0. (29) 

at 

By Gronwall’s inequality, we obtain that 

^Q{Qt,Pt,Zt) < e~'^^EQ{qo,po,zo), for all f > 0. 

We can re-write the above inequality using dehnition of {qt,pt,Zt) as follows 

EQ{{Q„Pt,Zt) - (Q,,Pt,%)) < e-^^EQ{{Qo,Po,Zo) - {Q,,Po,Zo)), 

for all t > 0. Now we optimize over {Qq, Pq, Zq) and {Qq, Pq, Zq) with re¬ 
spective laws po and pg to get 

EQ(((5o, -Po, Zq) — (Qg, P g, Zg)) = 1F2q(Po, Pq)^- (30) 

Then using the relation W 2 Q{pt,Pt)'^ < EQ((Qt, Pt, Z^) - {Q^,Pt,Zt)), we 
have 

^2Q{Pt,A) < VF2q(Po, Pg). (31) 

This completes the proof of the proposition. □ 

The following lemma will be used later to prove the existence and uniqueness 
of the stationary measure. 

Lemma 3. Lemma 7.3] Let (Al,dist) be a complete metric space and 
S{t) be a continuous semigroup on (Al,dist). Assume that there exists 0 < 
L{t) < 1 such that 

dist{S{t){x),S{t){y)) < L{t)dist{x,y), 

for all t > 0, and x,y in AA. Then there exists a unique stationary point 
Xoo e AA, i.e., S{t){xoo) = for all t > 0. 
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We are now ready to present the proof of Theorem [T] 

Proof of TheoremUl Let Q{q,p, z) be the positive dehnite quadratic form on 
obtained in Proposition [1] Since Q is equivalent to the canonical form 
kP + bP + bb there exists positive constants C and C such that 

< C"W2Q{pt,h) < C"e-^^W2Q{po,Po) < W2(po,Po), 

for alH > 0 and all solutions {pt)t>o and (j)t)t>o of ([I]) by Proposition [T] This 
proves the hrst assertion (jS]) of Theoremdl Since \/Q is a norm on R^'^, which 
is equivalent to the Euclidean norm \/\q\‘^ + bP + bP, according to [^, the 
space 1 ^ 22 ) is a complete metric space. The existence of a unique 

of stationary measure poo is then followed from the contraction property of 
Proposition [T] and Lemma O Moreover, taking {pt)t>o = Poo in (E]) we obtain 

dSD- □ 

3.2. Particle approximation of the main equation: proof of Theorem\^ 

We start with proving that the second moments of all solutions of ([T]) are 
hnite provided that the second moment of the initial data is hnite. 

Lemma 4. Under Assumption [U for all positive a, (3 and X there exists a 
positive constant po such that, if 0 < Ca + Cb < Po, then supi>o + 

bP + \^\‘^)dpt{q, P, z) is finite for all solutions {pt)t>o of O with initial data 
Po e 

The assertion of this lemma can be obtained using Theorem [H Indeed, 
let Poo be the unique stationary measures of ([I]) obtained in Theorem [TJ For 
all solutions {pt)t>o of ([I]) with initial data po G P 2 (R-^'^), from Theorem [H 
we have for all t > 0 

[ (bP + bP + \z\‘^)dpt{q,P,z) = W2ipt,5oY 

<2(W2(pi,Poo)' + M"2b0,Poo)') 

< 2(C"e-^*W2(po,Poo)' + W"2bo,Poo)') 

< 2 (C" hP2(po, PooY + W2{do, PooY) ■ 

However, to show that the hniteness of the second moments, and hence the 
property of propagation of chaos of the system ([7]), is independent of the long 
time behaviour of the solutions, we provide a direct proof using again the 
coupling technique similarly as in Theorem [H This proof also demonstrates 
further the usefulness of the coupling technique. 
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A direct proof of Lemma^ using coupling technigue. Let pt be a solution of 
([T]) with initial data po in Let ai,-- - ,05 be positive constants, 

which will be specihed later on. Dehne 

2(g,P, z) = ai\q\^ + a 2 |p|^ + + 2aiq ■ p + 2a5q ■ z + 2p- z. (32) 


Then 


d 

dt 




Q{q,p,z)pt{dqdpdz) = / Q{q,p, z)dtpt dqdpdz 


'R3d 


/R3d 


= 2asd + 2 


= 2asd+ 2 


Q{q,p, z) ( - divq(ppi) + divp[(/3g + A{q) + B * pt - \z)pt] 


/R3d 


/R3d 


+ div^[(Ap + az)pt] + dqdpdz 

p ■ {aiq + atip + Osz) - {fdq + A(g) + S * pi(g) - \z) ■ {a 2 P + a^q + z) 

— {Xp + az) ■ {a^z + a^q + p) pt{dqdpdz) 

— jSa^ |gp + (ai — {3a2 — Xa^) q ■ p + {—(3 + Aa 4 — 005 ) q ■ z 

— (A — 04) IpI^ + (as + Xa 2 — Xa^ — a)p ■ z — {aas — A) \z\‘^ 

- (A{q) + B * pt{q)) ■ {a 2 P + a^q + z) pt{dqdpdz). 

(33) 

using Young’s 


We now estimate the last two terms involving A and B in 
inequality and assumptions on A and B. 

For the A-terms, we get 


^\'^2F 






-I- 




< 2(7^(021^1 + 04^1 
A i7yi((2a4 + 02 + l)|q'P + 02^1 


2 , |„|2 


+ 



zf) + {(3af\q\‘^ + A|pp + aaf\z\ 
(34) 
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For the S-terms, we obtain 


2 / {a2P + a^q + z) ■ B * pt{q) pt{dqdpdz) 

=—2 {a2P + a^q + z) ■ B{q — q')pt{dq'dp'dz')pt{dq,dp,dz) 

J'RP'i 

= - [a2ip- p') + ai{q - q') + {z - z')] ■ B{q - q')pt{dqdp'dz')pt{dq, dp, dz) 

<Cb [a^\q - q'\ + a2\p - p'\ + \z - z'\]\q - q'\pt{dq'dp'dz')pt{dq,dp,dz) 
J'RP'l 


<Cb 


(*) 

< Cb 


?R6d L 


'R3d 


aM - q'? + y (k - + \p- p'?) 

+ - d'? + k - z'\^) pt{dq dp'dz')pt{dq, dp, dz) 

2ai\q\^ + a2(|g|^ + \p\^) + {\q\^ + |;2nl pt{dq, dp, dz) 


= Cb [{2a4 +a 2 + l)\q\‘^ + a 2 \p‘^\ + Izl"^] pt{dq,dp,dz), (35) 

where to obtain (*) we have used the following estimate 


/ \q-q'l"^ ptidq'dp'dz')pt{dq,dp,dz) = / {{ql"^ + {q'l"^ - 2q ■ q') pt{dq'dp'dz')pt{dq,dp,dz) 

/R6d jR6d 

2 


= 2/ \q\‘^pt{dq,dp,dz)-2( qpt{dq,dp,dz) 

jR3d V J^3d 


/R3d 

<2/ \q\^pt{dq,dp,dz), 

jRSd 


and similarly for other terms involving \p — p'l"^ and \z — z'l"^. 

To proceed, we hrst choose ai,..., as such that the terms involving inner 
products in fl55]) vanish, i.e.. 


{ ai - I3a2 - Aos = 0, 

—f3 + Aa4 — aos = 0, (36) 

Os T Aa2 — Aos — o = 0. 
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Substituting fl51|) and (I5S]) into fl5^ . and setting rj := Ca + Cb, we get 


d f 

37 / Q{q,p,z)pt{dqdpdz) 


dt 


'R3d 


2,(i^d 


Q<4 ^9 J- 

— + — H- 

/3 A aos 


|Al(0)p 


/R3d 


- (/3a4 - (204 + 02 + 1)??)) kP - (A - 204 - a2V) \p\‘^ 

— {aas — 2X — p)\z\‘^ pt{dqdpdz). (37) 


Similarly as in the proof of Theorem [T], we can prove the existence of a 
positive constant po, depending only on the parameters of the equation, such 
that for all 0 < Ca + Cb < Po, there exist Oi,..., 05 that satisfy fl5BD and 


{ / 3 fl 4 — (204 + 02 + l)p > 0, 

A — 2o4 — a2p > 0, (38) 

aos — 2X — p, 


and that Q{q,p,z) is a positive dehnite quadratic form on It follows 

from fl37D that 


d 

dt 


Q{q,p, z)pt{dqdpdz) < Ci - C 2 / (|gp + \p\^ + \z\^)pt{dqdpdz) 


'R3d 


< Cl - Ca 



Q{q,p,z)pt{dqdpdz), 


for some positive constants Ci, C 2 and C^. Applying the Gronwall’s lemma, 
we obtain that 

sup / Q{q,p, z)pt{dqdpdz) < 00 , 
t>0 J-R3d 

provided that initially Q(q,P, z)po(dqdpdz) < 00 . This is equivalent to 
the statement that 


sup 

t>o 


'RSd 


(kP + bP + jzj^)pt(dqdpdz) < oc, 


if initially po has hnite second moment. This completes the proof of the 
Lemma. □ 
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We are now in the position to prove Theorem [21 

Proof of Theorem [H We recall that for each 1 < "i < A^, the law fit of 
{Ql, PJ:, Zl) is the solution to Eq. ([T]) at time t with the initial datum 
Po and that the two processes {Q], PJ:, Zl)t>o and , Pt’^, Zf^)t>Q are 

driven by the same Brownian motions and initial datum. Dehne zl) := 

- {Ql, PJ, Zf), then {ql,pl, zl) satishes the following SDE 

'dqi = p\dt, 

dpi = -fdql dt - {A{QiP) - A{Ql)) dt-P E7=i(^(Qt’^ - Ql'"") - B * fitiQD) dt + Xzl dt, 

I 

dzl = —Xp\ dt — Xazl dt, 

A<li,pi,4) = (O.O.O). 

(39) 

The underlying idea of remaining of the proof will be similar as that of 
Theorem [1] and Lemma jH We will hnd a positive dehnite quadratic form TX 
on such that the SDE (15^ is dissipative with respect to V. Let ai,..., 05 
be positive constants to be chosen later on. We dehne 

V{q,p, z) := ai|g|^ + a2|p|^ + aal^l^ + 2aiq ■ p + 2a^q ■ z + 2p- z. (40) 

Let {ql,pl,zl) be a solution to the SDE fl39l) . We have 

jV{ql,pl,zf) 

= + osktT + 204?! • Pt + ‘^(^^ql ■ zl + 2p\ ■ zl) 

= 2aiql ■ Pt + 2^4blP + ^a^pl ■ zl 
+ ( 203 ^ + 2a^ql + 2p\) ■ {-Xp\ - azl) 

+ {2a,pl + 2a,ql + • ( - ^ql - {A{QJ^) - A{Ql)) - - - Qj'^) 

i=i 

- B » tLiQi)) + \zf) 

= (2ai — 2X05 — 2j3a2) ql ■ pI + (205 — 2Aa3 — 2 q! + 2Aa2) pi ■ zl + (2Aa4 — 2aa5 — 2/3) ql ■ zl 

- 2/304 biT - (2A - 204) bil^ - (2^03 - 2A) 

- {2a2pl + 204^/ + 2zl) ■ {A{Ql’^) — A{Ql)) 

1 ^ 

- + 2 « 49 ; + H) ■ WQf - «•") - B * MQi))- ( 41 ) 

i=i 
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Similarly as in fll9p we have 

— {2a2pl + 204^4 + 24)(^(Qt’’^) “ ^(Qt)) ^ ‘^<^2CA\<lt\\Pt\ + 204 ( 7 ^ 1 ^ 41 ^ + 2CA\ql\\zl\ 

< a^cA^+ bin+204^^194?+ caM^+ bin 

= (02 + 204 + l)( 7 ^|g 4 p + a 2 CA\pl\‘^ + CaIzII"^. 

(42) 


Substituting (14^ into (I4T|) . we obtain 

< (2ai — 2Aa5 — 2/3a2) b ' + (2^5 — 2Aa3 — 2a + 2\a2) p\- zl + (2Aa4 — 2005 — 2/5) q\ ■ zl 

— (2/5a4 — (02 + 204 + \)Ca) bn ~ ~ ~ 2a2(7yi) btl^ ~ (2Q;a3 — 2A — Ca) btl^ 

1 ^ 

- + 2«49; + 24) . {B{Q';” - Qf") - B * ,.,(«»)• («) 

i=i 

We first choose oi,..., 05 such that the cross terms in the right-hand side of 
fl43|l vanish, i.e., 

{ 2cl\ — 2Aci5 — 2ci2/5 = 0, 

2 X 0,2 — 2Aci3 -|- 2o^ — 2 q; = 0, (44) 

2Aa4 — 2aa5 — 2/5 = 0. 

With this choice, fl43|) reads 

jT^{qlPa4) 

< -(2^04 - (02 + 204 + 1)<7^) bir - (2A - 204 - 2o2(7a) bil" - (2ao3 - 2A - Ca) IbT 

1 ^ 

- ^ E(2«2pi + 204b + 24) ■ {B{Qf - Qf) - B * UQ\))- (45) 

i=i 

Note that this estimate holds for every z = 1,..., A/. By averaging fj45ll over 
i we have 

jB{ql,p\,zl) 

< -(2/5o4 - (02 + 204 1)Ca) btP “ (2A - 204 - 2o2(7a) bt P “ (2ao3 - 2A - Ca) btP 

2 ^ 

- E(“2P: + ad +4)- (B{Qd - dd) - B * M)). (46) 

*,i=i 
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Next we estimate the last term in 0461) . Using the following identity 

we now decompose the last term in fj46ll into six terms and estimate them as 
follows. 

1) By symmetry and the Lipschitz property of B, we have 

N 

- EK ■ WQi’' - Q’t’") - B(Qi - «))) 

N 

= ~Y1 EIW!'" - «:) ■ (-SWl" - «■") - - Q*))! 

*4=1 

1 " 

= -5 E *[((«< - ) - {Qi - 01)) ■ {B(Qf - 0?") - B{Qi - 0 *))] 

*4=1 


< E E - Qi'"') - (Qi 


N 


C 


*4=1 

N 


^ E E wl" -- («■ 


4,1V 


* 4=1 



2 


2 


N 1^2 

= CBY.m?-CB^\Y.<A 

ij=l 2=1 

N 

<CBEEi9:P = CBiv*Ei9;i*. 

*4=1 
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2) By assumption on B and Young’s inequality, we have 

N 

-Y,E\pi- (BlQi" - - B{Ql - Qi))] 

1 " 

= -5 E E[(p; -i4)- b(qi - «))! 

*J =1 


C 


N 


\pl-MQt''-Qt"") - (Q\ -Qi)\ 


i,j=^ 

N 


^ ^ Y E[i|p: -rir + Ikq;'" - «•") - (q: - «)i 

tj=i 


Cs 


N N 

]v"E|p,y - 1 EEpJf + -1 Y®« 

i=i i=i 


<%iv"E[i«,y + |p,y] 


3) Similarly we obtain 

N ^ r 

- E - «'") - -SW! - ^ + \4n 

*J =1 


4) We continue with the term 


N 


2 a4E 


Qt 


E {B(Qt - Qi) 


i=i 


B*p<ra)) 


< ]V/?a 4 Ek;!" + e| E (B{Qi - Qi) - B * MQ',))\ 

^ 2 

= NDatE\qf + Ee[-B(<?; -Qi)-B* MQ',) 

+ ^ E®[(-®W’< - «) - B» p.(q;)) ■ (B(« -Qi)-B* n(Qi)) 

^ 3^k 

(47) 
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We need to analyse further the last two terms in (H7|) . Because B is 
odd, we have -B(O) = 0. Hence for any y G R'^ it holds that \B{ii)\ = 
\B{y) — B{0)\ < Cslyl- It implies that 


E 


B{Ql-Q{)-B*y,{Qi) 


< 2 


E\BiQl-Ql)\‘^ + E\B*y,{Ql)\^ 


<2Cl 


<8Cl 


< M, 


ElQl-Q{l‘^+ I \q-q'\^yt{q,P,z)yt{q',p',z')dqdpdzdq'dp'dz' 






\q\‘^Pt{q,P, z) dqdpdz 


(48) 


where in the last inequality above we have used Lemma 0] to obtain 
M = 16(7^ {€' H 2 (po 5 Poo)^ + ^ 2(^05 Poo)^)- Certainly M depends only on 
the initial second moment and the coefficients of the equation and but not 
on t or N. 

The last term in flT7|l vanishes since for all j ^ k, we have 


E 


{B{Ql -Qi)-B* yM)) ■ {B{Ql - Qf) - B * ^QD) 

{E^.[BiQl - Qi) - B * y,iQl)]) ■ {EQ,[BiQl - Q^,) - B * y,iQl)]) 


- Egi 




(49) 


The equality above holds true because Ql and are independent and 
have the same law, which is the hrst marginal of pf 
Substituting flT 8 |) and flT9l) into fITTD . we get 


—2a4E 


N 

il ■ WQ\ - Qi) 

i=i 


BtMQi)) 


< iV/3o4E|<,;|2 + ^ M, 


from which, by summing over i, we obtain 

N 

-2a, Y, E Wt • {B{Q\ - Qi) - B * ptm) 


*j=i 


< N^lda,E\ql\^ + ^MN 

P 

= N^Pa,E\ql\^ + ^MN. 

P 

(50) 
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5 ) Similarly we obtain the following estimate 


N 


-2a2 Y, E [pl • {B{Qi - Qi) - B * pt{Ql)) 
i,j=^ 

6 ) Finally, for the last term we also get 


< N'^XE\pI\‘^ + ^MN. 

A 


N 


-2 5; E [j; . (B(Q‘ MQi)) 

*J = 1 


< N^aa3E\pl\^ + - MN. 

aa^ 


Set p := Ca + Cb- Bringing all terms together, it follows from 
d 


that 


dt 


z]) < -(/ 3 a 4 - (02 + 204 + l)h) \(l]? - (A - 204 - 2a2p) \p]\^ 

-(aa 3 - 2 A-r,)|z,y+(^+^ + ^)E. ( 61 ) 


Note the difference between the right hand side of the above inequality and 
that of (l 23 ll in Section Iddl The difference is due to the extra terms involving 
B{Q\ — Ql) — B*ptiQt) Cl'S shown in the above computations. We now choose 
04 = A. From (Hill we can express oi, 02 and 05 in terms of 03 and choose 03 
such that 

^ Aa4—/3 A^/4—/? 

“5 — a ~ a ’ 

02 = as + f - x> 


( 52 ) 


oi — Aos + /3 o 2 — Aos + /3 ^03 + y ~ j 

/ 3 o 4 — (02 + 204 + 1) 7 / = — (03 + y — ^ + A+ 1 ) 7 ]> 0 , 

A - 204 - 0277 = A - ^ j 7^ > 0, 

0)03 — 2A — 77 > 0, 


and that T){q,p, z) is a positive dehnite quadratic form on R^'^, where 

= (a5(A-y)+/3(o3+y)j l?l^+(«3+^-y) \p\‘^+a3\z\'^+^ q-p-2a5 q-z+2p-z. 

Similarly as in the proof of Theorem [Hand from (jhTll . we obtain the existence 
of a positive constant 770 depending only on a, / 5 , A, Ca and Cb such that for 
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all 0 < Ca + Cb, there exist oi,..., as such that V{q, p, z) is a positive dehnite 
quadratic form on and such that 

^EV(gl,pl zl) < -C'iE[|gip + \p]\^ + \zl\^] + (53) 

for alH > 0 and for positive constants Ci and (^ 2 , which depend only on the 
parameters of the equation and the initial second moment, but not on N. 
Since P is a positive quadratic on R^'^, the right-hand side of fl55D in turn is 
bounded by —C^iEV^ql^p], zj) + for some positive constant C^, so that 


, pI, zl) < -CsEViql, p \, zl) + ^. (54) 

From this inequality, we deduce that 

EF>(g 4 \pi\ zl) < for all t > 0, 

for some positive constant Using the fact that D is a positive dehnite 
quadratic on R^'^ again, we obtain the following bound 


eIIQ;" - Qil" + If, 


l,N 


pr+ iz 


l,N 





< 


c 

iV’ 


for all f > 0 and for some positive constant C depending only on the pa¬ 
rameters of the equation and the initial second moment, but not on N. This 
completes the proof of flTOl) . 

The assertion that the particle system ([7]) satishes the property of propa¬ 
gation of chaos is a direct consequence of ffTOj) . In deed, let (p,..., 4) be an 
arbitrary fc-uple of [1, A^] and (p be a 1-Lipschitz function on (R^*^)^. Then 


((R3d)fe J{R3d)fe 


ikX\ 






< 








, ? • • • 5 A^ 


,Xt 


< - x;X\ < 


ik-lX j^ikX^ 

Ck 


-E 


p{Xi-\...,X, 


+ ... 

ikX) 


i=l 
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where we have used the Lipschitz property of (p and ffTU]) . This estimate 
implies that converges to in the 1-Wasserstein distance, which 

is tested against 1-Lipschitz functions. Since the 1-Wasserstein distance is 
equivalent to narrow convergence plus convergence of the hrst moments, this 
concludes that p®!’"'’*'' converges to narrowly. 

Finally, the estimate flTT]) is followed straightforwardly thanks to fITU]) 
and definition of the Wasserstein distance m- We complete the proof of 
Theorem |2l □ 

References 

[1] L. Ambrosio, N. Gigli, and G. Savare. Gradient flows in metric spaces 
and in the space of probability measures. Lectures in Mathematics. ETH 
Zurich. Birkhauser, Basel, 2nd edition, 2008. 

[2] J. Bajars, J. Frank, and B. Leimkuhler. Stochastic-dynamical ther¬ 
mostats for constraints and stiff restraints. The European Physical Jour¬ 
nal Special Topics, 200(1):131-152, 2011. 

[3] J. Bajars, J. E. Frank, and B. J. Leimkuhler. Weakly coupled heat 
bath models for Gibbs-like invariant states in nonlinear wave equations. 
Nonlinearity, 26(7):1945-1973, 2013. 

[4] F. Bolley. Separability and completeness for the Wasserstein distance. 
In Seminaire de probabilites XLI, volume 1934 of Lecture Notes in Math., 
pages 371-377. Springer, Berlin, 2008. 

[5] F. Bolley, J. A. Gafiizo, and J. A. Garrillo. Stochastic mean-field limit: 
non-Lipschitz forces and swarming. Math. Models Methods Appl. Sci., 
21(11):2179-2210, 2011. 

[6] F. Bolley, A. Guillin, and F. Malrieu. Trend to equilibrium and particle 
approximation for a weakly selfconsistent Vlasov-Fokker-Planck equa¬ 
tion. M2AN Math. Model. Numer. Anal, 44(5):867-884, 2010. 

[7] J. A. Garrillo and G. Toscani. Gontractive probability metrics and 
asymptotic behavior of dissipative kinetic equations. Riv. Mat. Univ. 
Parma (1), 6:75-198, 2007. 

[8] D. A. Dawson and J. Gartner. Large deviations from the McKean-Vlasov 
limit for weakly interacting diffusions. Stochastics, 20(4):247-308, 1987. 


24 


[9] D. A. Dawson and J. Gartner. Large deviations, free energy functional 
and quasi-potential for a mean field model of interacting diffusions. 
Mem. Amer. Math. Soc., 78(398):iv-|-94, 1989. 

[10] M. H. Duong, M. A. Peletier, and J. Zimmer. GENERIC formalism 
of a Vlasov-Fokker-Planck equation and connection to large-deviation 
principles. Nonlinearity, 26(11):2951-2971, 2013. 

[11] J. Frank and G. A. Gottwald. The Langevin limit of the NosGHoover- 
Langevin thermostat. J. Stat. Phys., 143(4):715-724, 2011. 

[12] R. Kupferman. Fractional kinetics in kaczwanzig heat bath models. 
Journal of Statistical Physics, 114(1/2), 2004. 

[13] B. Leimkuhler, E. Noorizadeh, and F. Theil. A gentle stochastic thermo¬ 
stat for molecular dynamics. Journal of Statistical Physics, 135(2):261-- 
277, 2009. 

[14] S. Meleard. Asymptotic behaviour of some interacting particle systems; 
McKean-Vlasov and Boltzmann models. In Probabilistic models for non¬ 
linear partial differential equations (Montecatini Terme, 1995), volume 
1627 of Lecture Notes in Math., pages 42-95. Springer, Berlin, 1996. 

[15] M. Ottobre and G. A. Pavliotis. Asymptotic analysis for the generalized 
Langevin equation. Nonlinearity, 24(5):1629-1653, 2011. 

[16] A. A. Samoletov, C. P. Dettmann, and M. A.J. Chaplain. Ther¬ 
mostats for slow conhgurational modes. Journal of Statistical Physics, 
128(6):1321-1336, 2007. 

[17] A. Sznitman. Topics in propagation of chaos. In Paul-Louis Hennequin, 
editor, Ecole d’Et de Probabilits de Saint-Plour XIX 1989, volume 1464 
of Lecture Notes in Mathematics, pages 165-251. Springer Berlin Hei¬ 
delberg, 1991. 

[18] D. Talay. Stochastic Hamiltonian systems: exponential convergence to 
the invariant measure, and discretization by the implicit Euler scheme. 
Markov Process. Related Pields, 8(2): 163-198, 2002. Inhomogeneous 
random systems (Cergy-Pontoise, 2001). 


25 



[19] Cedric Villani. Topics in optimal transportation, volume 58 of Graduate 
Studies in Mathematics. American Mathematical Society, Providence, 
RI, 2003. 

[20] Cedric Villani. Hypocoercivity. Mem. Amer. Math. Soc., 
202(950);iv+141, 2009. 


26 



